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Introduction. 

The study of asymptotic (as independent variable tends to a singular 
point) properties of Painleve transcendents is one of the most important fields 
in modern theory of integrable nonlinear ODE's. The Painleve equations are 
known to be integrable in the sense of commutative matrix representation 
(Lax pairs). One has six matrix equations 

D,Lj - D,Aj + [Lj, Aj] = 0, J = 1, 2, 6, (1) 

where = d/dx; Lj = Lj{y,y',x,z), Aj = Aj{y,y',x,z) are 2*2 matri- 
ces that rationally depend on spectral parameter z, and the j-th Painleve 
equation y" — Pj{y,y',x) = is equivalent to (|^). The matrices Lj,Aj were 
written in paper 

The goal of this paper is to analyze asymptotic behavior of the sixth 
Painleve transcendent using the so-called Whitham method. The PVI case 
is tedious due to large amount of calculations, so it is easier to illustrate the 
basic ideas of the method (which are the same for all the six equations) on 
technically the simplest case of PI. 

The matrices Li and Ai look as follows: 

L - ( ^ ^] A - ( 2y + Az\ 

^'~[y-z J' ^'~[-x-y^ + 2yz-Az^ y' J' 

Introduce now new variable X and replace all the variables x explicitly en- 
tering formula (0) by X : Lj = Lj{y, y', X, z), Aj = Aj{y, y', X, z). For such 
matrices we have the following Lemma. 

Lemma 1. Let e be some real positive number. Then system 

eD.Lj - D^Aj + [Lj, Aj] = (3) 
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is equivalent to system 

D,X = e, y" = P^{y,y',X)=0. (4) 

Proof can be obtained via direct computation. So, for PI the system (^) 
has the form 

D^X = e, y"-3y^-X = 0. 

Calculations for the other Painleve equations are principally analogous and 
can be extracted from paper 0. 

Lemma 2. Solution of equation as e = and X = const can be 
represented by the following formula: 

l/o(x) = /,(r + J = 1,2,. ..,6, (5) 

where r = xU, U = U{a); fj are periodic functions which can be explicitly 
written out in terms of Weierstrass or Jacobi elliptic functions for any of the 
six Painleve equations. The vector a{X) consists of parameters that detrmine 
the elliptic function fj. $ is some phase shift. 

The proof uses the latter equation of system (H) where X is put to con- 
stant value. In the case of the first Painleve function fi is Weierstrass 
p— function: 

/i = 2p{x + $; g2, 93) ; 92 = -X, 93 = -Fi/4, (6) 

where Fi is some parameter, the formula (P) was first figured out in paper 
i- 

Now admit that number e ispositive and small. We look for solutions to 
equation in the form of formal series in parameter e : 

yix) = yo{x) + eyi{x) + (7) 

so that parameters determining the elliptic function yo = fj obey some special 
nonlinear ODE usually called Whitham equation or modulation equation. 
Thus, we look for the main term of series (0) in the form 

yo{T, X) = fj {e-'S{X) + <I>(X); a{X)) , DxS = U. 

Lemma 3. The Whitham equation can be written in the following form: 

Dx det Aj = a22DJii + auDj22 - ^12^^21 - 021^2^12, (8) 
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where Aj = (am„) , Lj = (Imn) , m,n = 1,2, the bar means averaging over 
period of the elhptic function (^. 

Proof. One can easily see that equation as e = indicates indepen- 
dence for spectral characteristics of matrix Aj of variable x. So the condition 

det Aj = holds. Formal introduction of variable X induces the change 
for differentiation rule: —>■ UDr + eDx, where parameter e is put to be 
small and positive. Further the condition (||) yields equation 

<m = e^^^n,m + [Lj, Aj]^,^ , 71,171 = 1,2. 

Substituting this into equality 

Dx det Aj = CiwCi22 ~l~ '-^22'^ll — 0,^2^21 — 0,2iO,i2, 

we change the differentiating rule and obtain the following: 

{UDr + eDx) det Aj = e ia22DJu + aiiDj22 - aiiD.ki - 021^^/12) +0 (e^) , 

Now average, i.e. integrate over the period (in "fast" variable r). The aver- 
aging kills complete derivatives in r which gives the claim. 

Corollary 1. There exists unique coefficient of the polynomial det Aj with 
non-trivial dynamics in X in force of the modulation equation. Denote this 
coefficient Fj. Thus the Whitham system can be written as unique ODE on 
Fr 

The corollary can be verified via direct calculations for all the six equa- 
tions. For PI we have the following: 

detAi = 16z^ + 4Xz - Fi, 

where Fi = (y')^ — 2y^ — 2yX. The modulation equation (^) takes the form 

Dx det Ai = 4:Z + 2y 

and can be rewritten as DxFi = —2y. Taking into account the solution (|^), 
we obtain: 

DxFi = -2r]/u = 2ei + 2(63 - ei)E/K, 
where E = E{k), K = K{k) are complete elliptic integrals: 



K=l . E= f\r-^'''dz. k^='-^ 



1- z^il-k^z^) Jo ^ l-z^ ei-es 
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ei,2,3 are roots of Weierstrass polynomial Rsit) = it^ — g2t — gs; g2 = 
-X, gs = -Fi/4. 

Corollary 2. The simplest way for obtaining the elliptic ansatz fj is solving 
equations 

Fj = consti, X = const2. (9) 

Lemma 4. The elliptic ansatz @ forms the main term in series in small 
parameter e (0) for solution to system (^. 

To prove this one should see that perturbation of solution to system @) 
with e = runs continuously while e obtains small non-zero value. The 
appropriate elementary calculations are illustrated here on the simplest ex- 
ample of PI. So, for e > system (|]) is y" = Sy"^ + Xq + ex, where Xq is 
constant. Via simple manipulations this can be reduced to condition 

2dx = , , = + 0(e) 

\'2y^ + Xq + const 

which means that the main term of the series (|^) is the function fi (see (^) 
on condition that x does not belong to small neighborhoods of singularities 
of the elliptic function /i. 

Now we can prove the following theorem. 

Theorem 1. The function ?/o determined by formulas (H) and (H) forms the 
main term of asymptotic series for solution of appropriate Painleve equation 
as |x I tends to infinity. 

Proof. The scale transformation x —>■ ex leads to change D^X = e 
DxX = 1 in formula (^. Therefore the expansion (|^) in small parameter e 
turns to series in negative powers of large variable x. 

2. PVI and the Whitham method. 

The sixth (and the most common) Painleve equation 

. 1 (i , 1 , ^) U/f JU^ ,^) y', (10) 

^V^ y ~ ^ y ~ ^ J \^ X — 1 y — X J 

y{y — l){y — x) / _|_^^_|_ x—1 

x^{x — iy \ y'^ {y — ly {y ~ ) '' 
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where the Greek letters denote free parameters, can be obtained as the com- 
patibihty condition of the following linear system of equations: 

D,Y^Ae{z,x)Y{z,x), D^Y ^ Le{z,x)Y{z,x), (11) 

where 

A(zx)^( "^^(^'"^^ "^^(^'"^^ ]=^ + ^ + JL 
' \a2i{z,x) a22{z,x) J z z-1 z - x' 

^' - ( !' «.) ) . - 0. Ia(.. .) ^ -A'^^. (12) 



Put 



fci 

fc2 



012(2;) = 



ki + k2 = - {00 + 01 + e^) , ki-k2 = 9oc, 
luqUq uJiUi UxUx k{z - y) 



z z — 1 z — x z{z — l){z — x)' 



u^au{y)^ \ — H , (13) 

y 1/ - 1 y-x 

u = -a22(y) 



X 



y y-1 y-x 
Then uq + Ui + Ux = k2, coqUo + ooiUi + oOxUx = 0, 

Uq + 9o Ui + 9i Ux + Ox , , , , , 

1 1 = 0, [X+ljCUoUo+XUJiUi+UJxUx — K, XCUqUo — K{x)y, 

UJq UJi LUx 

which are solved as 

_ ky _ k{y - 1) _ k{y - x) 

'^0 — ) '^1 — ~1 ) 



XUq^ [x — l)Ui x{x — l)Ux^ 

y Q 

where 

5i = y{y -l){y- x)u^ + [e^{y - x) + xOxiy - 1) - 2^2(1/ - l){y - x)] u+ 



5 



my - X - 1) - /cs (^1 + xO^) , 

where 

Si = y{y-l){y-x)u''+[{ei + e^){y - x) + xd,{y - 1) - 2k2{y - l){y - x)] u+ 

kl{y-x)-k2{ei+xe^)-kik2, (14) 
^ ^ 1/ - a; ^ 
a;(x - l)6'oo °°' 

where 

5oo = -x) + x{9^ + ^oo)(y - 1) - 2k2{y - l){y - x)] u+ 

kl{y -1) - k2 {01 + xOnc) - xkik2. 
The compatibihty condition for ( pTj) imphes 

L_h^J^y^\. (15) 

a;(a;-l) V V-^ y-xj 

Thus ?/ satisfies PVI with the parameters 

Now we apply ideas described in the previous paragraph to asymptotic 
analysis of the sixth Painleve transcendent. First calculate determinant for 
the matrix Aq. Using formulas ([l3|) , ([l^ one obtains: 

011(2;) = —R~^{z)S^ 

where 

S = kiz"^ + z[x (mo + Ml + ^0 + ^1) +uo + eo + u^ + e^] + o (2°) , 

0'22{.z) = -R^^{z) \k2Z^ - z\{x + 1)Uq + XUi +Ux] + (-2°)} , 

where R{t) = t{t—l){t—x)] O (z^) means powers of z of order not higher than 
j. The entries ai2, (221 yield terms of lower order in z, so they can be ignored 
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while computing the two higher terms of polynomial det Aq. Therefore we 
have the following: 

detAe = R^\z)S, 

where 

S = kik2z'^— 

[ki {uo{x + 1) + UiX + Ux) - k2 {x {uo + Ui + 60 + 61) + Uq + 9o + Ux + Ox)] + 

0(z'' 



Setting 

det Aq = R-\z) [kik2Z^ + F^z^ + O [z^)] (16) 
we get the coefficient Fq that determines the Whitham dynamics: 

Fe = {ki - k2) {ui + xux) - x {2kik2 + 9x) - 2kik2 - fc2^i. (17) 

The current goal is to extract the constraint on elliptic function from 
condition (|T7|). To do this we use (IH): 



6*00 {ui + xUx) = -R{y)u^+ 

u [{y - l)(y - x) {2k2 + ^^oo) -ei{y-x)- xOxiy - 1)] + 
kik2ix + I - y) + k2 iei+x9x), 
which via (^) and (|T^ turns to the following: 

, \ x'^ (x — 1)^ , .-2 

Ooo [ui + XUx) = — ^^^^^ [y ) + 

^y'x{x-l) |b + -J-^ [(y - l){y - x) {2k2 + ^00) - 0,{y -x)- xOxiy - 1)]| - 

^R{y)B^ + [xOxiy -l)-{y-l){y- x) (2^2 + ^oo) +Oi{y- x)] + 
kik2ix - y + 1) + k2 {ei+x9x), 

where 

„ ^0 , 01 Ox + i 

B = \ H . 

y y-l y-x 
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Now substitute this into ( ]T7| ) and obtain final constraint on genus one Rie- 
mann surface (y', y) and appropriate elliptic uniformization (consider x and 
Fq parameters): 



x\x - If (y'f - 2y'x{x - l)yiy - I) + y^ \l - (k^ - hY 
2y^ [(fci + ^2) C - 1 + 2x6, (1 - k2) + 2Fq] - 

y's+ 

2yx [2kik2ix + 1) + 2x6, (1 - ^2) + 2^6 - 60C] - x^6l = 0, 
where C = (x + 1) {ki + k2) + x6, + 6*1, 



+ 



(18) 



S = C^-l-2x6o {ki + k2)+4:kik2 (x^ + x + l)+4a;(a;+l)(l-A;2)^.+4(x+l)F6 

To start the Whitham asymptotic analysis we need also the modulation 
equation in addition to ansatz (0). It can be found in the following way. 
First replace variables x to X in formula (|l^) and differentiate it in X : 



Dx det Ae = ^r^^ {'2kik2 + DxF^) + O (z^) . 



(19) 



On the other hand we have condition (j|) which is to be studied now. Thus 
we have the following: 



D^L^{z,x) 



A. 



whence obtain: 



0'llDzl22 



{z-xy 



DJ 



z''22 



{z-xy 



a22Dzl 
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R{z){z-Xy 

Ux + 6, 
R{z){z-Xf 



z^ki + 0{z 
Z^k2 + 0{z) 



substitute into (||) and get: 



Dx det Aq 



z^ [u, {ki - k2) - k26,] + 0{z) 
R{z){z-Xf 



(20) 
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where bar means averaging. Comparison of formulas ([T^) and (^0]) yields the 
modulation equation: 

DxFe = {h - k2) - he^ - 2kik2. (21) 

One can as well rewrite ( ^I]) in the initial coordinates y, X. To do this use 

ih - h) = j^l'^^f , (22) 

where 



S = R{y)u^ + u[eiiy-X)+X {6, + 9^) - 1) - 2k2{y ~l){y- X)] + 
kl{y -l)-k2 {Oi + Xe^) - Xk,k2. 



Now substitute (|22|) into (0), again utilize (0), (0), (P^, simplify and 
finally obtain the modulation equation: 

D.F, 4 (fa - k) D,y + - f^^>^-J^ + (23) 



X{X - 1 



-5+ 



2(x - 1; 



[^0 {h - h) + 2X (2A:iA;2 + d^) + 2^2 (^i + ^2 + ^i) + 2Ff, 
-k29^ - 2kik2, 



where 



S = -{k2- k^) [X {k2 - A;i - 6^) + 6, + 6, + 1] - 

X {2kik2 + e^) - k2 {ki + k2 + 9i) - Fq. 

Here y denotes the mean for elliptic function y specified by equation 
where Fg and X (instead of x) are considered as parameters. 

3. Partial solutions for the modulation equation and PVI. 



Analysis of the system (|T8|), (p3D in generic form is cumbersome, moreover 
there is a question of the phase shift $ within the elliptic ansatz (^. This is 
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why we start with the simplest solutions that correspond to strongly degen- 
erate surface (piSl) . While trying to find partial solution to ([T8|) - (|23|) among 
elementary functions one can note asymptotic homogeneity of formula (|I^) 
for large x. Denote y = and rewrite ([T8| ) in variables x^' and ^. The 
discriminant for this polynomial looks as follows: 



D = i^ {k2 - hf - 2^' [{k2 + k^f + 29, (1 - k2) + 2^60;-^ + O [x'^) 
[(^2 + kif + Akik2 + AO, (1 - k2) + AFqx-^ + O (x-^ 



(24) 



2^ 2F6a;-2 + 0(x 



Oix 



Seeking condition for strong degeneracy of the Riemann curve ( 
out that the polynomial (p4D tends to have two double roots if 
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and Fq = -2kik2X + o{X), X ^ oo, 



one find 



(25) 



D = e{^-ink2-k^f + 0{X- 

In this case four branch points of the Riemann surface ( [T^ ) asymptotically 
coinside pairwise (that is what we call double or strong degeneracy). Substi- 
tuting condition ( P3| ) into (|THp one easily obtains the appropriate asymptotics 
for solution: ^ = 1 + o(l) and, therefore, 



y = X + o{x), X oo, (26) 
where o{x) denotes terms that grow not faster than log x. One can also easily 



verify that solution (^SQ , (^q) suits the modulation equation (|23|) . 
So we have the following theorem. 

Theorem 2. In the case of 9, = {S = 1/2) the sixth Painleve equation 



has a solution with asymptotics (|26|)p|. 

To prove this one should note in addition to mentioned above that strong 
degeneracy of the elliptic ansatz (|T^) transforms the phase shift $ in formula 
(^ into a shift in variable x which can be found via simple iterative procedure 
computing terms of the series 



^Such a solution for PVI under a 
H was found in paper B 



(2/i - l)/2, /? = 7 = 0, (5 = 1/2 and half-integer 
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